Clover: An unbiased method for prioritizing differentially expressed genes using a data-driven approach.

Identifying key genes from a list of differentially expressed genes (DEGs) is a critical step in transcriptome analysis. However, current methods, including Gene Ontology analysis and manual annotation, essentially rely on existing knowledge, which is highly biased depending on the extent of the literature. As a result, understudied genes, some of which may be associated with important molecular mechanisms, are often ignored or remain obscure. To address this problem, we propose Clover, a data-driven scoring method to specifically highlight understudied genes. Clover aims to prioritize genes associated with important molecular mechanisms by integrating three metrics: the likelihood of appearing in the DEG list, tissue specificity, and number of publications. We applied Clover to Alzheimer's disease data and confirmed that it successfully detected known associated genes. Moreover, Clover effectively prioritized understudied but potentially druggable genes. Overall, our method offers a novel approach to gene characterization and has the potential to expand our understanding of gene functions. Clover is an open-source software written in Python3 and available on GitHub at https://github.com/G708/Clover.

Morc1 reestablishes H3K9me3 heterochromatin on piRNA-targeted transposons in gonocytes.

To maintain fertility, male mice re-repress transposable elements (TEs) that were de-silenced in the early gonocytes before their differentiation into spermatogonia. However, the mechanism of TE silencing re-establishment remains unknown. Here, we found that the DNA-binding protein Morc1, in cooperation with the methyltransferase SetDB1, deposits the repressive histone mark H3K9me3 on a large fraction of activated TEs, leading to heterochromatin. Morc1 also triggers DNA methylation, but TEs targeted by Morc1-driven DNA methylation only slightly overlapped with those repressed by Morc1/SetDB1-dependent heterochromatin formation, suggesting that Morc1 silences TEs in two different manners. In contrast, TEs regulated by Morc1 and Miwi2, the nuclear PIWI-family protein, almost overlapped. Miwi2 binds to PIWI-interacting RNAs (piRNAs) that base-pair with TE mRNAs via sequence complementarity, while Morc1 DNA binding is not sequence specific, suggesting that Miwi2 selects its targets, and then, Morc1 acts to repress them with cofactors. A high-ordered mechanism of TE repression in gonocytes has been identified.

Protocol for identifying differentially expressed genes using the RumBall RNA-seq analysis platform.

Here, we present a protocol for the identification of differentially expressed genes through RNA sequencing analysis. Starting with FASTQ files from public datasets, this protocol leverages RumBall within a self-contained Docker system. We describe the steps for software setup, obtaining data, read mapping, sample normalization, statistical modeling, and gene ontology enrichment. We then detail procedures for interpreting results with plots and tables. RumBall internally utilizes popular tools, ensuring a comprehensive understanding of the analysis process.

Modeling embryo-endometrial interface recapitulating human embryo implantation.

The initiation of human pregnancy is marked by the implantation of an embryo into the uterine environment; however, the underlying mechanisms remain largely elusive. To address this knowledge gap, we developed hormone-responsive endometrial organoids (EMO), termed apical-out (AO)-EMO, which emulate the in vivo architecture of endometrial tissue. The AO-EMO comprise an exposed apical epithelium surface, dense stromal cells, and a self-formed endothelial network. When cocultured with human embryonic stem cell-derived blastoids, the three-dimensional feto-maternal assembloid system recapitulates critical implantation stages, including apposition, adhesion, and invasion. Endometrial epithelial cells were subsequently disrupted by syncytial cells, which invade and fuse with endometrial stromal cells. We validated this fusion of syncytiotrophoblasts and stromal cells using human blastocysts. Our model provides a foundation for investigating embryo implantation and feto-maternal interactions, offering valuable insights for advancing reproductive medicine.

Loop-extruding Smc5/6 organizes transcription-induced positive DNA supercoils.

The structural maintenance of chromosomes (SMC) protein complexes-cohesin, condensin, and the Smc5/6 complex (Smc5/6)-are essential for chromosome function. At the molecular level, these complexes fold DNA by loop extrusion. Accordingly, cohesin creates chromosome loops in interphase, and condensin compacts mitotic chromosomes. However, the role of Smc5/6's recently discovered DNA loop extrusion activity is unknown. Here, we uncover that Smc5/6 associates with transcription-induced positively supercoiled DNA at cohesin-dependent loop boundaries on budding yeast (Saccharomyces cerevisiae) chromosomes. Mechanistically, single-molecule imaging reveals that dimers of Smc5/6 specifically recognize the tip of positively supercoiled DNA plectonemes and efficiently initiate loop extrusion to gather the supercoiled DNA into a large plectonemic loop. Finally, Hi-C analysis shows that Smc5/6 links chromosomal regions containing transcription-induced positive supercoiling in cis. Altogether, our findings indicate that Smc5/6 controls the three-dimensional organization of chromosomes by recognizing and initiating loop extrusion on positively supercoiled DNA.

Cell cycle-dependent gene networks for cell proliferation activated by nuclear CK2α complexes.

Nuclear expression of protein kinase CK2α is reportedly elevated in human carcinomas, but mechanisms underlying its variable localization in cells are poorly understood. This study demonstrates a functional connection between nuclear CK2 and gene expression in relation to cell proliferation. Growth stimulation of quiescent human normal fibroblasts and phospho-proteomic analysis identified a pool of CK2α that is highly phosphorylated at serine 7. Phosphorylated CK2α translocates into the nucleus, and this phosphorylation appears essential for nuclear localization and catalytic activity. Protein signatures associated with nuclear CK2 complexes reveal enrichment of apparently unique transcription factors and chromatin remodelers during progression through the G1 phase of the cell cycle. Chromatin immunoprecipitation-sequencing profiling demonstrated recruitment of CK2α to active gene loci, more abundantly in late G1 phase than in early G1, notably at transcriptional start sites of core histone genes, growth stimulus-associated genes, and ribosomal RNAs. Our findings reveal that nuclear CK2α complexes may be essential to facilitate progression of the cell cycle, by activating histone genes and triggering ribosomal biogenesis, specified in association with nuclear and nucleolar transcriptional regulators.

Churros: a Docker-based pipeline for large-scale epigenomic analysis.

The epigenome, which reflects the modifications on chromatin or DNA sequences, provides crucial insight into gene expression regulation and cellular activity. With the continuous accumulation of epigenomic datasets such as chromatin immunoprecipitation followed by sequencing (ChIP-seq) data, there is a great demand for a streamlined pipeline to consistently process them, especially for large-dataset comparisons involving hundreds of samples. Here, we present Churros, an end-to-end epigenomic analysis pipeline that is environmentally independent and optimized for handling large-scale data. We successfully demonstrated the effectiveness of Churros by analyzing large-scale ChIP-seq datasets with the hg38 or Telomere-to-Telomere (T2T) human reference genome. We found that applying T2T to the typical analysis workflow has important impacts on read mapping, quality checks, and peak calling. We also introduced a useful feature to study context-specific epigenomic landscapes. Churros will contribute a comprehensive and unified resource for analyzing large-scale epigenomic data.

FXYD3 functionally demarcates an ancestral breast cancer stem cell subpopulation with features of drug-tolerant persisters.

The heterogeneity of cancer stem cells (CSCs) within tumors presents a challenge in therapeutic targeting. To decipher the cellular plasticity that fuels phenotypic heterogeneity, we undertook single-cell transcriptomics analysis in triple-negative breast cancer (TNBC) to identify subpopulations in CSCs. We found a subpopulation of CSCs with ancestral features that is marked by FXYD domain-containing ion transport regulator 3 (FXYD3), a component of the Na+/K+ pump. Accordingly, FXYD3+ CSCs evolve and proliferate, while displaying traits of alveolar progenitors that are normally induced during pregnancy. Clinically, FXYD3+ CSCs were persistent during neoadjuvant chemotherapy, hence linking them to drug-tolerant persisters (DTPs) and identifying them as crucial therapeutic targets. Importantly, FXYD3+ CSCs were sensitive to senolytic Na+/K+ pump inhibitors, such as cardiac glycosides. Together, our data indicate that FXYD3+ CSCs with ancestral features are drivers of plasticity and chemoresistance in TNBC. Targeting the Na+/K+ pump could be an effective strategy to eliminate CSCs with ancestral and DTP features that could improve TNBC prognosis.

Context-dependent perturbations in chromatin folding and the transcriptome by cohesin and related factors.

Cohesin regulates gene expression through context-specific chromatin folding mechanisms such as enhancer-promoter looping and topologically associating domain (TAD) formation by cooperating with factors such as cohesin loaders and the insulation factor CTCF. We developed a computational workflow to explore how three-dimensional (3D) structure and gene expression are regulated collectively or individually by cohesin and related factors. The main component is CustardPy, by which multi-omics datasets are compared systematically. To validate our methodology, we generated 3D genome, transcriptome, and epigenome data before and after depletion of cohesin and related factors and compared the effects of depletion. We observed diverse effects on the 3D genome and transcriptome, and gene expression changes were correlated with the splitting of TADs caused by cohesin loss. We also observed variations in long-range interactions across TADs, which correlated with their epigenomic states. These computational tools and datasets will be valuable for 3D genome and epigenome studies.

Nutrient-driven dedifferentiation of enteroendocrine cells promotes adaptive intestinal growth in Drosophila.

Post-developmental organ resizing improves organismal fitness under constantly changing nutrient environments. Although stem cell abundance is a fundamental determinant of adaptive resizing, our understanding of its underlying mechanisms remains primarily limited to the regulation of stem cell division. Here, we demonstrate that nutrient fluctuation induces dedifferentiation in the Drosophila adult midgut to drive adaptive intestinal growth. From lineage tracing and single-cell RNA sequencing, we identify a subpopulation of enteroendocrine (EE) cells that convert into functional intestinal stem cells (ISCs) in response to dietary glucose and amino acids by activating the JAK-STAT pathway. Genetic ablation of EE-derived ISCs severely impairs ISC expansion and midgut growth despite the retention of resident ISCs, and in silico modeling further indicates that EE dedifferentiation enables an efficient increase in the midgut cell number while maintaining epithelial cell composition. Our findings identify a physiologically induced dedifferentiation that ensures ISC expansion during adaptive organ growth in concert with nutrient conditions.

Phase-separated nuclear bodies of nucleoporin fusions promote condensation of MLL1/CRM1 and rearrangement of 3D genome structure.

NUP98 and NUP214 form chimeric fusion proteins that assemble into phase-separated nuclear bodies containing CRM1, a nuclear export receptor. However, these nuclear bodies' function in controlling gene expression remains elusive. Here, we demonstrate that the nuclear bodies of NUP98::HOXA9 and SET::NUP214 promote the condensation of mixed lineage leukemia 1 (MLL1), a histone methyltransferase essential for the maintenance of HOX gene expression. These nuclear bodies are robustly associated with MLL1/CRM1 and co-localized on chromatin. Furthermore, whole-genome chromatin-conformation capture analysis reveals that NUP98::HOXA9 induces a drastic alteration in high-order genome structure at target regions concomitant with the generation of chromatin loops and/or rearrangement of topologically associating domains in a phase-separation-dependent manner. Collectively, these results show that the phase-separated nuclear bodies of nucleoporin fusion proteins can enhance the activation of target genes by promoting the condensation of MLL1/CRM1 and rearrangement of the 3D genome structure.

Comprehensive multiomics analyses reveal pervasive involvement of aberrant cohesin binding in transcriptional and chromosomal disorder of cancer cells.

Chromatin organization, whose malfunction causes various diseases including cancer, is fundamentally controlled by cohesin. While cancer cells have been found with mutated or misexpressed cohesin genes, there is no comprehensive survey about the presence and role of abnormal cohesin binding in cancer cells. Here, we systematically identified ∼1% of cohesin-binding sites (701-2,633) as cancer-aberrant binding sites of cohesin (CASs). We integrated CASs with large-scale transcriptomics, epigenomics, 3D genomics, and clinical information. CASs represent tissue-specific epigenomic signatures enriched for cancer-dysregulated genes with functional and clinical significance. CASs exhibited alterations in chromatin compartments, loops within topologically associated domains, and cis-regulatory elements, indicating that CASs induce dysregulated genes through misguided chromatin structure. Cohesin depletion data suggested that cohesin binding at CASs actively regulates cancer-dysregulated genes. Overall, our comprehensive investigation suggests that aberrant cohesin binding is an essential epigenomic signature responsible for dysregulated chromatin structure and transcription in cancer cells.

Mesenchymal loss of p53 alters stem cell capacity and models human soft tissue sarcoma traits.

Soft tissue sarcomas (STSs) are a heterogeneous group of tumors that originate from mesenchymal cells. p53 is frequently mutated in human STS. In this study, we found that the loss of p53 in mesenchymal stem cells (MSCs) mainly causes adult undifferentiated soft tissue sarcoma (USTS). MSCs lacking p53 show changes in stem cell properties, including differentiation, cell cycle progression, and metabolism. The transcriptomic changes and genetic mutations in murine p53-deficient USTS mimic those seen in human STS. Furthermore, single-cell RNA sequencing revealed that MSCs undergo transcriptomic alterations with aging-a risk factor for certain types of USTS-and that p53 signaling decreases simultaneously. Moreover, we found that human STS can be transcriptomically classified into six clusters with different prognoses, different from the current histopathological classification. This study paves the way for understanding MSC-mediated tumorigenesis and provides an efficient mouse model for sarcoma studies.

Transient Methionine Deprivation Triggers Histone Modification and Potentiates Differentiation of Induced Pluripotent Stem Cells.

Human induced pluripotent stem cells (iPSCs) require high levels of methionine (Met). Met deprivation results in a rapid decrease in intracellular S-adenosyl-methionine (SAM), poising human iPSCs for differentiation and leading to the apoptosis of undifferentiated cells. Met deprivation triggers rapid metabolic changes, including SAM, followed by reversible epigenetic modifications. Here, we show that short-term Met deprivation impairs the pluripotency network through epigenetic modification in a 3D suspension culture. The trimethylation of lysine 4 on histone H3 (H3K4me3) was drastically affected compared with other histone modifications. Short-term Met deprivation specifically affects the transcription start site (TSS) region of genes, such as those involved in the transforming growth factor ß pathway and cholesterol biosynthetic process, besides key pluripotent genes such as NANOG and POU5F1. The expression levels of these genes decreased, correlating with the loss of H3K4me3 marks. Upon differentiation, Met deprivation triggers the upregulation of various lineage-specific genes, including key definitive endoderm genes, such as GATA6. Upon differentiation, loss of H3K27me3 occurs in many endodermal genes, switching from a bivalent to a monovalent (H3K4me3) state. In conclusion, Met metabolism maintains the pluripotent network with histone marks, and their loss potentiates differentiation.

CohesinDB: a comprehensive database for decoding cohesin-related epigenomes, 3D genomes and transcriptomes in human cells.

Cohesin is a multifunctional protein responsible for transcriptional regulation and chromatin organization. Cohesin binds to chromatin at tens of thousands of distinct sites in a conserved or tissue-specific manner, whereas the function of cohesin varies greatly depending on the epigenetic properties of specific chromatin loci. Cohesin also extensively mediates cis-regulatory modules (CRMs) and chromatin loops. Even though next-generation sequencing technologies have provided a wealth of information on different aspects of cohesin, the integration and exploration of the resultant massive cohesin datasets are not straightforward. Here, we present CohesinDB (https://cohesindb.iqb.u-tokyo.ac.jp), a comprehensive multiomics cohesin database in human cells. CohesinDB includes 2043 epigenomics, transcriptomics and 3D genomics datasets from 530 studies involving 176 cell types. By integrating these large-scale data, CohesinDB summarizes three types of 'cohesin objects': 751 590 cohesin binding sites, 957 868 cohesin-related chromatin loops and 2 229 500 cohesin-related CRMs. Each cohesin object is annotated with locus, cell type, classification, function, 3D genomics and cis-regulatory information. CohesinDB features a user-friendly interface for browsing, searching, analyzing, visualizing and downloading the desired information. CohesinDB contributes a valuable resource for all researchers studying cohesin, epigenomics, transcriptional regulation and chromatin organization.

HSF1 phosphorylation establishes an active chromatin state via the TRRAP-TIP60 complex and promotes tumorigenesis.

Transcriptional regulation by RNA polymerase II is associated with changes in chromatin structure. Activated and promoter-bound heat shock transcription factor 1 (HSF1) recruits transcriptional co-activators, including histone-modifying enzymes; however, the mechanisms underlying chromatin opening remain unclear. Here, we demonstrate that HSF1 recruits the TRRAP-TIP60 acetyltransferase complex in HSP72 promoter during heat shock in a manner dependent on phosphorylation of HSF1-S419. TRIM33, a bromodomain-containing ubiquitin ligase, is then recruited to the promoter by interactions with HSF1 and a TIP60-mediated acetylation mark, and cooperates with the related factor TRIM24 for mono-ubiquitination of histone H2B on K120. These changes in histone modifications are triggered by phosphorylation of HSF1-S419 via PLK1, and stabilize the HSF1-transcription complex in HSP72 promoter. Furthermore, HSF1-S419 phosphorylation is constitutively enhanced in and promotes proliferation of melanoma cells. Our results provide mechanisms for HSF1 phosphorylation-dependent establishment of an active chromatin status, which is important for tumorigenesis.

CRISPR/Cas9 Screening for Identification of Genes Required for the Growth of Ovarian Clear Cell Carcinoma Cells.

Epithelial ovarian cancer is classified into four major histological subtypes: serous, clear cell, endometrioid and mucinous. Ovarian clear cell carcinoma (OCCC) responds poorly to conventional chemotherapies and shows poor prognosis. Thus, there is a need to develop new drugs for the treatment of OCCC. In this study, we performed CRISPR/Cas9 screens against OCCC cell lines and identified candidate genes important for their proliferation. We found that quite different genes are required for the growth of ARID1A and PIK3CA mutant and wild-type OCCC cell lines, respectively. Furthermore, we found that the epigenetic regulator KDM2A and the translation regulator PAIP1 may play important roles in the growth of ARID1A and PIK3CA mutant, but not wild-type, OCCC cells. The results of our CRISPR/Cas9 screening may be useful in elucidating the molecular mechanism of OCCC tumorigenesis and in developing OCCC-targeted drugs.

Large-scale multi-omics analysis suggests specific roles for intragenic cohesin in transcriptional regulation.

Cohesin, an essential protein complex for chromosome segregation, regulates transcription through a variety of mechanisms. It is not a trivial task to assign diverse cohesin functions. Moreover, the context-specific roles of cohesin-mediated interactions, especially on intragenic regions, have not been thoroughly investigated. Here we perform a comprehensive characterization of cohesin binding sites in several human cell types. We integrate epigenomic, transcriptomic and chromatin interaction data to explore the context-specific functions of intragenic cohesin related to gene activation. We identify a specific subset of cohesin binding sites, decreased intragenic cohesin sites (DICs), which are negatively correlated with transcriptional regulation. A subgroup of DICs is enriched with enhancer markers and RNA polymerase II, while the others are more correlated to chromatin architecture. DICs are observed in various cell types, including cells from patients with cohesinopathy. We also implement machine learning to our data and identified genomic features for isolating DICs from all cohesin sites. These results suggest a previously unidentified function of cohesin on intragenic regions for transcriptional regulation.

Cohesin-dependent chromosome loop extrusion is limited by transcription and stalled replication forks.

Genome function depends on regulated chromosome folding, and loop extrusion by the protein complex cohesin is essential for this multilayered organization. The chromosomal positioning of cohesin is controlled by transcription, and the complex also localizes to stalled replication forks. However, the role of transcription and replication in chromosome looping remains unclear. Here, we show that reduction of chromosome-bound RNA polymerase weakens normal cohesin loop extrusion boundaries, allowing cohesin to form new long-range chromosome cis interactions. Stress response genes induced by transcription inhibition are also shown to act as new loop extrusion boundaries. Furthermore, cohesin loop extrusion during early S phase is jointly controlled by transcription and replication units. Together, the results reveal that replication and transcription machineries are chromosome-folding regulators that block the progression of loop-extruding cohesin, opening for new perspectives on cohesin's roles in genome function and stability.

PHOSPHATE exporter XPR1/SLC53A1 is required for the tumorigenicity of epithelial ovarian cancer.

Ovarian cancer is the fifth most common cause of cancer-related death in women. Ovarian clear cell carcinoma (OCCC) is a chemotherapy-resistant epithelial ovarian cancer with poor prognosis. As a basis for the development of therapeutic agents that could improve the prognosis of OCCC, we performed a screen for proteins critical for the tumorigenicity of OCCC using the CRISPR/Cas9 system. Here we show that knockdown of the phosphate exporter XPR1/SLC53A1 induces the growth arrest and apoptosis of OCCC cells in vitro. Moreover, we show that knockdown of XPR1/SLC53A1 inhibits the proliferation of OCCC cells xenografted into immunocompromised mice. These results suggest that XPR1/SLC53A1 plays a critical role in the tumorigenesis of OCCC cells. We speculate that XPR1/SLC53A1 might be a promising molecular target for the therapeutic treatment of OCCC.